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ACCOMPLISHMENTS 


An oral presentation was made on July 23 to Mr. F. W. Steinle 
and personnel of the Experimental Investigations Branch, NASA Ames. 

At that time, a review of technical progress was given. 

As discussed in detail within the accompanying report, a correla- 
tion of Preston-tube data with theoretical skin-friction coefficient 
has been achieved for the subsonic, compressible laminar boundary 
layers on the AEDC Cone. The recommended correlation has been 
developed using data from nineteen different wind-tunnel conditions 
and has an rms error in skin-friction coefficient of less than 5%. 

The STAN-5 computer program for boundary layer calculations is 
not sensitive to changes in cone pitch or yaw angles. Thus, if the 
effect of such angles on the correlation is to be studied, a more 
sophisticated analysis of three-dimensional, viscous flow will be 
needed, e.g., McRae, et al. [Ref. 1]. 

The simplified model for calculating the magnitude of Preston- 
tube pressures {as a function of boundary-layer profile, local static 
pressure, and probe geometry and position with respect to the wall) 
does not appear to be a fruitful approach. Thus a more rigorous 
analysis will be necessary if this type of sensitivity study is to 
be physically meaningful. 

An approach for developing a correlation for the subsonic, 
turbulent boundary layer and transitional region has been selected. 
Skin friction and velocity profiles, at the beginning of the 
turbulent boundary layer, can be estimated by using the correlation 
of Allen [Ref. 2] in conjunction with the Preston-tube data and 
the Wu and Lock and STAN-5 computer programs. Once the distribu- 
tion of turbulent skin-friction and boundary layer profiles are 
available, a correlation between Preston-tube data and theoretical 
skin friction can be developed using the same techniques employed 
for the laminar boundary layer. Skin friction within the transi- 
tion zone can be easily approximated by employing the empirical 
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interim ttency function of Dharwan and Narasimha [Ref. 3], Although 
this i ntermi ttency functioi is based on flat-plate measurements, the 
use of actual Preston-tube measurements to specify the extent of the 
transition zone will result in a very good approximation for the dis- 
tribution of C f through the transition zone. 

In the case of laminar boundary layers, there is no need to 
employ the more sophisticated program of Wilcox and Rubesin. However, 
this program may still be useful in checking the STAN-5 results for 
compressible, non-adiabatic turbulent boundary layers. This analysis 
and_option will be relegated to future work. 

The supersonic wind-tunnel data cannot be successfully analyzed 
without a calibration of P ref as a function of Preston-tube position, 
M «> and Re ft< The corresponding calibrations for the flight experi- 
ments could conceivably be utilized, but this analysis will also be 
relegated to future work. 
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REMAINING TASKS TO BE ACCOMPLISHED 
UNDER THIS GRANT 

The efiects of changes in nose bluntness on pressure distribu- 
tion along the AEDC cone will be investigated. 

Subsonic Preston-tube data will be used to study and compare the 
onset and extent of boundary layer transition for the corresponding 
flight and wind-tunnel flow conditions. 

U^e the flight data to develop a correlation for subsonic laminar 
boundary layers, with and without heat transfer, and compare the 
results with the corresponding correlation of the wind-tunnel data. 

The pressure distribution, measured during flight, will be used to 

calculate the flow, rather than the theoretical pressures of Wu and 
Lock. 
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abstract 

The laminar ooundary layer on a 10-degree ccne ir a tran- 
sonic wina tunnel Is studied. Tne Inviscic flow and boundary 
layer development are simulated Dy computer programs. The 
effects ot pitch and yaw -ingles on_ the boundary lay^r are 
examined. 

Rr estori- tube data, taken cn the Arnold engineering 
evfclopment Center CAEDC) boundary-Layer-Transition Cone In 
the NASA Ames 11-it Transonic wind Tunnel, has been used tc 
develop a correlation which relates tne measurements to 
theoretical values ol laminar skin friction. The recommended 
correlation is cased on a compressible torn ol the classical 
1 aw-of-the-wail . 

The computer codes successfully simulate the laminar 
bouncary layer tor near-zero pitch ana yew angles. However, 
In cases of significant pitch and/or yaw angles, the flow is 
three-dimensional and the boundary layer computer code useo 
iiuri: cannot provide a satisfactory model. 

The skin-friction correlation is thought to be valid 
for body geometries other than cores. It accounts lor varia- 
ble pi op «c ty anu heat transfer effects. Tne rms deviatior 
between tneoretical skin-friction coefficients ano the cor- 
responding correlation values is < 5 t. Thus, 


as perhaps 


mi.,ht o« expected, this is * hatter correlation tor 
compressible iaainar flows than has been reported for com- 
pressible, turbulent layers. The new correlation can be 
employed in tr an so nlc-wlno- tunnel tests to relate preston- 
tube surveys along models to distributions oi laminar, 
sKin-f ricti on coefficient. 
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nomenclature 


focal skin-friction coefficient 

Specific heat at constant pressure,* 0.24 Utu/iboi°H 
tor air 

Preston- tuo e pressure coef i Iclent, =( P n . -p )/(0.5pU 2 ) 

Ft w ' e e 

-C-har acter is tic dimension ot the probe, in. 
equivalent circular diameter ot the probe, in. 
hlasius velocity ratio,* u/n 

e 

Gain factor lor the pressure transducer, psi/m. 
Converter factor,* 32.174 !bm-ft/lbf-s 


Enthalpy, Utu/lbm 
Pressure head, in. 

Mechanical to thermal energy 
77H.2 ibt-tt/Ptu 

Non-oiiuo/tsional normal distance,* 
Cone axial length, in. 

Mach numoer 
Pressure, psi 
Pranoti number 


conversion factor,* 


2y/n 
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Dynamic pressure, pst 

Recovery factor, or racial distance. In. 
Gas constant,* 53 . 35 1M -f t/lbin°R tor air 
heynulds nuniuet based on n ar.d U e D/ v w 

Freestream unit Reynolds number,* U.^/v^ 


x 


Ro l.»n*jth Reynolds numbur/- II x/ t , 

x e e 

Re„ Moinentum-tnickness Reynolds number/* U 6/v 
0 e w 

T Temperature/ 0 !? 

T* Reference temperature/ 0 !? 

u Longitudinal velocity Inside boundary layer/ ft/s 

Up^ Mean velocity across probe face/ ft/s 
u* Sheai veloc ity/ = Vr ™//T 

u + Normalized velocity for wail-laW/= u/u* 

U Velocity outside boundary layer/ ft/s 

x Distance along cone surface/ ft 

X Body force per unit volume/ ibt/ft 

x* D ins ens ion less independent variable/ fign. (8.13) 

y Distance normal to cone surface/ ft 

y* Dimensionless dependent variable/ Kqn. (8.8b) 

y hall Reynolds number/* y u*/v w 

SuUwCimid 


aw adiabatic wall 

B HldSius solution 

c tor cone flow 

e at ecgt of boundary layer 

eff effective 


eq equivalent 

E external or outer 

I Internal or inner 

Pt Pros ton-tub e 

ref reterence 

s shorted 


JJIUGINMS * 
of. VOOU OU.M.nY 


t total 

w at wail or cone surface 

W tor wedge flow 

o wind-oif 

°* treestraaro condition 
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P 
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8 * 
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V 
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v' 
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* 

to 


Angle of attack/ deg. 

Effective angle of attack/ deg. 

Yaw angle/ — fte_g. , or pressuce gradient paraueter 
katlo of specific heats, = 1.4 for air 
Cone semi-vertex angle, deg., or bounoary layer 
thickness, ft 

Displacement thickness of boundary layer, it 
Deflection or increment 
Aziroutn angle, deg. 

blasius non-dimensional normal distance, = y^U e /2 xi/' 
Momentum thickness ot bounoary layer, ft 
Molecular viscosity, Ibf-s/ft 
Kinematic viscosity, ft?/s 

Kinematic viscosity evaluated at the reference temp- 

o 

erature, ft /s 
Density, lorn/ ft 
Shear stress, psf 

Angle between cone axis and resolved yaw vector, deg. 

Stream function, Ibm/s 

Normalized stream function, Ejn.(6.1) 
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CHAPTFK I 


INTRODUCTION 

The overall oojective ot this research is a better under- 
standing of boundary lay.er transition/ as reflected in the 
capabiLl_ty to relate transition on models in tr ansoniX— ulnd 
tunnels to the corresponding free-f light conditions* The 
particular objective of the work reported herein is to 
develop a correlation which relates Prestont-.tube measure- 
ments witnin a laminar boundary layer on a cone tc the cor- 
responding theoretical values ot skin friction. 

Pr eston- tub e measurements along the surface of a sharp 
l-O-deyr ee cone were obtained in the NASA Ames 11-Ft Tran- 
sonic wind Tunnel 03. The minimum and maximum pressure 
locations/ obtained during a survey along the length of the 
cone/ were interpreted as the onset and end of transition/ 
respective! y . 

The oouncary layer on the slender cone was simulated 
via the STAN^ computer code f2l which is an extended version 
ot latankar and Spalding's noundary layer program C3j. The 
invisctd tlow was calculated with Viu and Lock's computer 
program CO/ and the results were used as boundary condi- 
tions along the outer eaqe ol the bounaary layer. Subrout- 
ines were added to this program so that arbitrary combina- 
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tl0P " 01 l ' ttrn drd y ' lh an « ,es h « input, and tne pressure 
'list rihu tAon along the ray corresponding to the Preston-tube 

survey is always generated, in adcltion, a subroutine was 

added to the Wu and Lock program to calculate the initial 
Protlies needed for STANS, 

the cone is assumeo to be stationary, smooth ano 
snarp-nosed. 71, e probe is assumed to be stable, in contact 
with -the cone surlace, and lie totally inside the boundary 
layer, trie flu, is assumed to be ail-symmetrlc, adiabatic, 
compressible and althout body forces. The llo. outside the 
boundary layer ls assumed to be invlscid and is calculated 
based on the cone geometry, ,.e., viscous Interaction Is 

Ignored. The study was restricted to laminar hcuncary layers 
on the cone at subsonic speeds. 

The eftect on the invlscid nos of yarn and pitch ancles 
less than the cone seml-verte* anyle is easily calculated 
“itt. tae „u and Lock program. Hovevec, tne STAN5 program Is 
a t.o-1, mansions, boundary layer cede and uas found to be 
relatively insensitive to cnangos in these angles. 

* least-squares curvo-Uttiny program C101 .as used to 
•nr.ve at e Simple correlation between skin friction ano 

Pres tun- tube measurements lol the laminar, subsonic boundary 
1 dyfir. 


CHAPTER IT 


OBJFCTIVES 

The first objective ot this study was to calculate the best 
possible Initial profiles, which are required to begin num- 
erical o ouncsry- 1 aye r calculations, so that boundary-layer 
predictions would be uniformly accurate. In an earlier study 
by huprikar C51, it was found that different starting pro- 
files resulteo in differences in the computed shear stress 
near the tip ot the cone. An example of this is shown in 
Figure 1. 

The second objective was to extend the functions of Wu 
and Lock's program 4.41, which calculates tne inviscid pres- 
sure distribution on sharp cones at transonic Mach numbers, 
so as to automate calculation of tne pressure distribution 
along a ray correpondiny to the Preston-tube survey for 
norwero pitch ( ot ) and yaw (0) angles. This information 
then provides tne lnvincio boundary conditions for calcula- 
tion of tne boundary layer with STAN5. The thiro objective 
was to obtain a correlation for skin-friction coefficient or 
wall-shear stress in terms of tne Preston-tube pressure mea- 
surements, so that the Preston tube can be used as a skm- 
frlction measuring device. 

ORIGINAL PACK IS! 

OF POOR QUALITY) 
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The present research focuses on the NASA Ames wind tun- 
nel data taken within laminar boundary layers on tne AEDC 
Transition Cone at subsonic speeds. 


ORIGINAL PAGE IS 
OE POOR QUALITY, 




CHAPTER III 


EXPERIMENTAL DATA 


3.1 General Background 


The "ensure., nts utm.ee ln_ms research mere obtained i„ 
the NASA IX-Pt Transonic bind Tunnel at Noflet Field, Call- 
transcntc elhd tunnel Is an experleental taclllty 
intended to simulate the Uo« o,er scaled, aerodynamic-test 
models that uould be similar to lull-scale eenlcles during 

free-tllgnt through the atmosphere at Kach nuabers tro. 
^^pro^lmately 0.5 to 1.5. 

in transonic floe the duterence be.u.en the treestrea. 
velocity ano the speed of sound Is small compared to the 
magnitude ot either, and the chants in these parameters are 
Ot comparable magn.tude. This ,s contrasted to subsonic 
flo-e -here the velocity is loser than the sonic speed anc 
-nere Causes each number are primarily due to changes ,r 
creestreae velocity at essentially constant some speeds, 
and to supersonic Uos shore the magnitude o, the treestrea. 
velocity is substantially larger than the local sonic spaac 
slth changes in Hach hummer securing through variations ol 

froth parameters. f n t h e transomr m .. v 

me transonic Mach number range, not 

only oo cofrpresatollity eifects become 1 


mportant, compared 
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1 

to lower subsonic ttacn numbers where tne tlow is 
incompressible/ but also the tlow at near-sonic speeds is 
complex because of the mixed type or flow which may exist 

1 

with local supersonic flow tields contained in subsonic flow ( 

regions or local subsonic flow lielos embedded in supersonic 
flow regions* That is why the cone shape is used as a model 
for boundary-layer-transition research; since it will not 
have local shocks along the conical surface. At high sub- 
sonic speeds/ a shock may be generated near the base of the 
cone owing to flow expansion at the rear of the conical sur- 
face and a subsequent recompress ion in the wake. At super- 
sonic speeds/ tne core will/ of course/ also generate a oow 
shock/ but a shock does not occur on the surface throughout 
tne subsonic Mach number range. j 

it is worth mentioning that tne ventilated/ test-sec- 
tion wails cf a transonic wind tunnel introduce accustic ano : 

streamline disturbances into the test-section flow -which I 

means that the wind tunnel flow does not correspond exa ct l y 

i 

to transonic/ free-tlight conditions C63. No satisfactory \ 

\ 

1 

method has yet been derived to correct for all of tne wall j 

effects C71/ although this is an area of active research 

on. 

i.2 Apparatus end Measurements 

I 

Tne experimental data were obtained irom a Pitot probe that 1 

was traversed longitudinally along the surface of a 5-degree 
hall-angle cone. The cross-section of the opening of the 


H 


probe is shown In Flyura 2 loj. The opening has on oval 
shat>e witn the 3 raall dimension nornai to the cone surfaca. 
The outer heignt ot the probe face is 0.0097"/ while the 
centerline of the opening is 0.004S3" above the cone sur- 
face. A schematic of the experimental model and instrumenta- 
tion is shown In Figure 3. 

The total pressure/ as sensed by the Fitot probe#, was 
measured by a differential pressure transducer. The refer- 
ence pressure for the transducer was taken from tne static 
holes on a flow-angularity probe mounted underneath the 
cone. 

The output from the pressure transducer/ AH/ was 
recorded/ during constant wind tunnel conditions/ as a func- 
tion of x cn a plotter. Snorted output of the transducer/ 
for the same wind tunnel conditions/ was also plotted on the 
s d me plot. The output ot the transducer/ when tne tunnel was 
off and tne transducer was shortec/ was also plotted. This 
output should theoretically be zero. This deflection is 
called * wind-off* deflection. 

Using this information/ tho total pressure Ppt' mea- 
sured by the Pitot probe can be deduced by using the rela- 
t ion 

p Pt = P ref + G(AH + AH s + AH o> l3,U 

here F ref la the reference static pressure whicn is consid- 
ered to be egu.U to the freestream static pressure E9J. AH 


tho del lection oi the plotter corresponding to tne magni- 
tude o£ AP aa aensea by the differential pressure trans- 
ducer. The deflection from the shorted output Is AH S , and 
AH q is the vand-oft deflection. G is the gain factor of the 
Plotter/ and its value is 0.2bl5 psi/in. This value was det- 
ermined froa the calibration of the plotter E53. 

Twenty-one cases were chosen for detailed analysis. 
These are ail the available, subsonic-wind-tunnel cases with 
with hoar-zero Liow angles. The tabulated data for these 
runs is s nown in Table I. The freestream Mach number, unit 

Reynolds number and dynamic pressure are given oy M , R e 

oo It 

,lmi q » ' respectively, while « and /? are the angles of 
attack and yaw, respectively. 
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TABLE 1 
CASES STUDIED 


RUN NO , 

14 

00 

Re ft x 10 " b 

9 * 

o ° 


15.231 

0.95 

4 

693 

- 0.048 

0.018 

19.289 

0.8 

4 

617 

- 0.003 

- 0.022 

21 .318 

0.7 

4 

548 

- 0.006 

- 0.025 

23 . 346 

0.6 

4 

477 

- 0.001 

- 0.025 

25.376 

0.5 

4 

404 

- 0.005 

- 0.025 

27.411 

0.4 

4 

403 

- 0.004 

- 0.026 

29.440 

0.3 

4 

230 

- 0.006 

- 0.026 

39.545 

0.4 

2.5 

396 

0.023 

0.021 

40.547 

0.6 

5 

586 

0.021 

0.021 

41.548 

0.7 

5 

680 

0.018 

0.021 

42.549 

0.8 

5 

761 

0.013 

0.021 

43.550 

0.9 

5 

842 

0.010 

0.021 

44.551 

0.95 

5 

873 

0.008 

0.021 

56.631 

0.9 

3 

492 

0.062 

0.006 

57.632 

0.8 

3 

453 

0.066 

0.006 

58.633 

0.7 

3 

408 

0.071 

0.006 

59.634 

0.6 

3 

357 

0.075 

0.006 

60.635 

0.5 

3 

302 

0.068 

0.007 

61.636 

0.4 

3 

246 

0.070 

0.007 

70.726 

0.7 

4 

538 

0.036 

0.023 

72.748 

0.8 

4 

605 

0.030 

0.023 
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CALCULATION PROCEDURE: 

The calculation procedure consists of the following steps 
for each case studied : 

1. The given freestream parameters ( M w /Re^/q^, )/ as well 
as the tlow angles ( a / fi ), are fed into the extended Wu 
and Lock program. This program is described in the next 

chapter and is listed In Appendix F. The output is two- 

fold: 

a. The invlscio velocity distribution along the cone/ and 

b. The initial protiles of velocity and stagnation 
enthalpy at a distance very close to the tip ot the cone. 

?• These results — axe-then input to the STAN5 program. A 
brief description of STAN5 is presented in Chapter VI. The 

output trotr that program has detailed irtormaticn on the 

boundary-layer properties along the ray o£ the cone which 
corresponds to the Preston- tube survey. 

when correlation ot skin friction was pursued/ two more 
steps were followed : 

3. Experimental hr es ton- tube pressure measurements were cal- 
culated from NASA/Anes 11 TWT Preston-tube data C 1 3 using 
" gua t ion (3*1). 
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14 

4. These experimental pressures/ together with some other 
parameters* calculated by STAN5/ Here ted into a curve-fit- 
ting program £101 to obtain the required correlation* 

Figure 4 is a flow chart that summarizes the calculation 
procedure described above. 


‘See Chapter VIII tor oetails 


M »* %,* R e ft , a, 0 



!• Effective center of probe 
2. Effect of a, 0 
. . . etc. 


figure 4. Flow Chart tor tha Analyst 







CHAPTER V 


LX TEN OL'D VU AND LOCK COMPUTER PROGRAM 
5.1 Introduction 

Wu and Lock C4J developed a computer program to calculate 
the inviscic transonic flow fielc over a sharp-edge smooth 
cone surface. The program appears to give accurate results 
191 when compared with experimental observations. The pro- 
gran/ however/ nancies only yaw angles less than the cone 
semi-vertex angle. The modified program presented herein 
calculates the following aduitional Information 5 

1. The mvlscld velocity arid pressure 
distribution along a ray of the cone that 
corresponds to the Preston-tube survey for 
arbitrary combinations of a and p , 

2 . Tne effect of yaw angles on the inviscic flow 
field/ jnd 

3. The velocity and total enthalpy profiles at a 
user-speci f led initial station. 

The listing oi the extended program can he found ir Appendix 
F. output for an example run (Case 25.376) is also lrcluded. 

5.? The Original Program 

The main program reads in M w / a* fit and the core semi- 
vertex angle 8 and and calculates the inv lsc ic-t low pressure 


i* 
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distributor along the cone surface. This is used ay the 
pressure at the edge ot the boundary layer. The theory ano 
equations used are described In the Mu and Lock report C43. 
The program prints out, along the cone length, the local 
Mach number M g and the ratio P w /P w . 

5.3 Subroutine ANGLES 

Tills subprogram -uses the angle-of-attack a anc the yaw 
angle ft to calculate the effective pitch angle 5 and the 
azimuthal position of the probe e . This subroutine uti- 
lizes the equations derived by Dunn et al CllJ. These equa- 
tions are presentee m Append* A. the pcjah.e position is 
considered to oe always at the top of the cone and. In 
accordance with W u and Lock's notation, e =C.O always cor- 
responds to the leeward side of the cone. The calculated 
angles 15 and e) are then used in the main program to cal- 
culate the irvtscid pressure distribution along the top of 
the cone. 


5. 4 Subroutine DIST 

Tills subroutine reads In the freestream dynamic pressure 
(UfNF), unit Reynolds number (RLFT) and Mach numoer (MINE). 
It then uses tnese values to calculate the freestream pro- 
perties CPINF,T1NK,RH0INF,MUINF) as well as the tctal temp- 
erature and pressure <TTOT,ptot). The equations used are the 
equation of state tor a perfect gas (air), Sutherland's 
equation of viscosity and the isentropic relations Lii), The 


lb 


details ot tho calculations are described iri Appencix B. 

Next/ tnc subroutine uses the local Mann rumbers at 
stations along the cone surface, which are calculated In the 
main program, to calculate the local temperatures and veloc- 
ities using lsentroplc relations. These velocities are then 
used as the outer boundary conditions for calculation ot the 
boundary layer using STANS. 

5.5 Subroutine tNITIA 

This last subprogram calculates the velocity and stagnation 
enthalpy profiles across the boundary layer at a specifies 
initial location. It calculates the average static tempera- 
ture and viscosity across the bouncary layer anc uses tnei 
to modify the fiat-plate Blasius solution so as to apply to 
the cone problem. Tne details ate presented In Appendix C. 

5.o Checking Hu anc Lock Calculations 

As a check on the reliability of our version cf Hu and 
Lock’s program, the inviscid flow was calculated for a 
lu-gpyree cone at a 2-oegree pitch angle ano compared with 
tnose in wu and Lock’s report C4J. The following observa- 
tions were made : 

<j, Static pressures on the wlrdwaid side ot the cone art 
larger than tnose on the leewaro sice. 

b. increasing a increases the static pressure on the 
windward sloe and decreases it on the leeward side. 


c. The slope of the pressure clstribution is essentially 


*ne s units 

on butu 

the 

windward 

and 

I eeward 

sides oi the 

cone 

( except 

near the 

tip 

and the 

rear 

ends ol 

the core). 

Thes 

cnecks a 

re shown 

in 

Figure 5* 
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Figure i>. Calculated Pressure Coefficient on Cone Sur- 
face for Various Pitch Angles 



CHAPTER VI 


STAN5 COMPUTKK PftOGRA w 

Hased oil tbe work ot Patankac and Spalding £33/ the STAN5 
code mas developed by Crawfore anc Kays £23 as an implicit/ 
finite-difference/ f orward-marching Integration procedure 
which may fee used for computer simulation of bouncary layers 
with transition. The program solves simultaneously equations 
for conservation of mass/ momentum/ stagnation enthalpy and 
up to five mass transfer equations. 

The program uses eitner two-oimensionai plar.ar or axi- 
symnietnc type of coordinates so that it is possible tc 
solve for a large variety ot rlows by simple manipulation of 
variables. Tills is accomplished by replacing the y-coorui- 
nate with the stream runction 1 1/ < The u-velccity component 
is ae fined by 


and the momentum and ereryy equations become 


pu If + pu 


3 __ 

dtp 


r p u ij 


9u 

eff 


= -g 


dP 
c dx 


+ 9 C X ’ 


,na pu |St + u|j 


r 2 p « 

Kr eff 


Iff pH ' 


[ p eff / 

1 J- \ 

k J \ 

Pr eff / 


2 3 

r p u w 


(-ft] 


n 


Tl 

T ii t stream uinct ion t p is then normal lied by using trip trnns- 
f ormation 

$ 

w a (6,1) 

where and ip j ar e the stream iunctton values on the boundary 

suriaces or boundary conditions. 

A micro-integral method is usee to obtain implicit fin- 
ite-dif f erence equations/ which model the partial differen- 
tial equations and nay be used in a downstream/ forward- 
marchiny solution scheme. The program solves laminar and 
turbulent boundary layers, boundary-layer transition is 
based on the rpotrentum-thiCKPess Reynolds number criterion/ 
wnich is deiinej a priori. The way the transition Reynolds 
number (RETkAN) is specified is as follows ; 

A very large value Is assigned to RKTRAN, e.g. 10000/ so 
that the program is ensured to run wholly laminar. From the 
experimental data sheets obtained from NASA ClJ/ the loca- 
tion oi tne minimum pressure is considered to be the onset 
of transition. At this location the corresponding value of 
Re fl in 5TAN5 output is then considered to be the correct 
RliTRAH. 

However/ since we ate presently concentrating only on 
rhe laminar boundary layer/ a large value ot reTKAN was 
always assigned in the input to STANb and no re-run was 


f i e c e s s a r y 


'n 


OUii’i input pa fame tors and ••ilagw" ace retired# a 
detailed description of which can no found in the STAN5 
report L3J. The edge velocity distribution and initial pro- 
files for the velocity and total enthalpy across the boun- 
dary layer are required input tor STANS. They are prepared 
by the extended Wu ana Lock program. (See Chapter V). 

The output ot the program gives/ at every incremental 
X/ all the boundary- layer properties of interest/ e.g./ 
u(y)/U+(y)/y+(y),C f / S/ 6*/ Q, R e<? /P w /T(y)/ T t (y),P t (y), ... 
etc. This in tor ma 1 1 on can then be used for theoretical ana- 
lysis ot the boundary layer. 


f. 
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CHAPTER yll 

PtFECT OF PL Ok ANGLES ON BOUNDARY-LAYER 
CALCULATIONS 

As Mentioned bet-ore in Chapter V/ tno angles of the frees- 
tream flow will affect the boundary-layer flow** One ot the 
objectives cf this research was to investigate the capabil- 
ity ot trie available computer programs (fcu and Lock*s ana 
T A N r > ) to handle pitch ana yaw angles that are a significant 
fraction of and to obtain some conclusions regarding the 
analytical tools needed to analyze such cases* 

The original Viu and Lock program was modified to calcu- 
late the eflectlve yaw angle for arbitrary combinations ot 
yaw and pitch aryles. Tne equations derived by Dunn et al 
r 11 J ware used in subroutine ANGLES to calculate the azirauth 
angle ot the L' r o h> « ^ as discussed in Chapter V. It was found 
tnat the extended ku and Lock program works well with all 
the cases studied* 

One case was studied in some detail/ viz./ Case 40. 
which Mas the following data : 

M,„ " 0-6, R« ft 5 x 10 6 , q o » 586 psi , 

« = B = 0 .Oi!l °. 


‘Transition is affected when u / is chanyec by + b%. 
(Reterencu 9). 





Tut;, case was pick no up ii :i a .‘itdrt becaust: of ft:; relatively 
'•mall Macu number would allow neylection oi noise etfuct 
Tne output ot the program (edge velocities and 
Initial profiles) was input to STAN*. it was found that the 
results ot STANb for this cast were exactly the same as the 
cast of zero flow angles* This was not unexpected since (a ) 
and (0) are wary small in this case. 

Then, the same case was repeated but with larger 
angles, viz., a *2.0, /3=2.0 oeyrees, which places the probe 
135° from the windward element, and a=-P.O, /3=2.0 degrees 

which places the probe 45° iron; the wineward element* The 
results of these two runs, together with the ongiral run, 
are snown In f'tgure f. k is defined as 2y/D. The plottec 
results agree witn the observation that the pressures on the 
wino war c sloe ere greater tnan those on tne leewarc side and 
that a l et n - inci do nc e flow lies In between these. However, 
by comparing the values of wall shear stress, at the probe 
azimuth angle/ and the boundary-layer thicknesses 8, §*, 

0, the effect of a and /3 is negligible as sho*r in table 
If. It was, therefore, decided to confine the present stage 
ol research tu the cases of very small How angles* A possi- 
ble reason for STANDS Insensitivity is Its assumption ol 
ax i symmetry while real liow with large pi ten ara/or yaw 
angles will fuv<i siyriil .leant cress iluw/ thus forming a 
tnr ef-dimens iona 1 flow. 


, psf 


2800 


270C 


2600 


2500 


M oo=0.6, Ref t = 5x10 6 

• Experimental 
D STAN5, e = 135° 

° STAN5, a = /3=0.021 
v STAN5, e = 45° 



}k =1.2 


table II 

SENSITIVITY OF STAN5 COMPUTATIONS TO 
CHANGES IN FLOW ANGLES 3 


= 6 = 0.021 


0 [c f /2 


0.1 842 615 153 56 

0.2 593 811 216 78 

0.3 483 999 265 95 1482 11020 

0.41417 I 1157 305 110 j 41 7 I 1182 

0.5| 374 I 1310 341 123 | 372 | 1322 

0.6| 341 | 1449 374 134 | 3401 1425 


0 C f /2 6 


153 56 843 574 

216 78 593 853 

264 95 484 1025 

305 110 419 1188 

341 122 375 1330 

374 134'-- 


6* 


153 

56 

217 

78 

265 

96 

306 

no 

342 

123 

374 

134 


All numbers, except x, are multiplied by 10 6 . The case 
ana yted is 40.547. The aaimuthal angle, c, corresponds 
to the probe located at the top of the cine corresp ° n<) ' 
























CHAPTER VIII 


COPhFLATION OK SKIN FkICTIOn 


0*1 Theoretical Background 

h iraension*i] analysis £171 has led tc a wall-law ol the term 

u + = f(y + ) 

+ + (8.1) 

wnere u = u/u*, y = u*y/v anQ 

u-lonyitudinal velocity, u* = shear velocity = VTTp 
U *° Us * irconpressiDle hernculli equation to relate 
u to the procure difference between a Preston tube rest- 
ing cn the wall and local static pressure, the law-ot-the- 
waii can then be written in the following form. 

u Pt ■ V'2(P Pt - P w )/p' -/Tp 7?' 

U " e Iurth,!r “»«• th*t the ef lective center ot the 
IU0he y y sff , 1„ j t | ts hil 1 1 -nc lcjh t, K eff , 2, ef f /DM . 0 , 

then 


“ ’(V") = g (faTJZ 


(8.4) 


°* fruit ip ly tne riuii-eratoi -«nj denominator on the lett by 
appropriate factor. In oroex to obtain the same ytoupmy 
01 mn,s fll-Poai in the function y. 

*n n 2 /e 2 / / t — ' 
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1*4 J 4Li£\ lfl.5) 
4p v" \4p v / 

Tnis last relation provides a convenient way for determining 
the skin friction since tne shear stress is new uniquely 
related to the difference in pressure head measured with a 
Prpston- tube static-hole combination. For a Preston tube ol 
qlven geometry/ the function F can theoretically oe deter- 
mined from pipe flow experiments where the skin friction can 
ho deduced from measurements of pressure drop. In view of 
the t ic t that the wall-laws for pipe and bouncary layer 
flows are identical C173/ the calibration is expected to 
hole also in oounoary layer flows. 

Equation (y.b) is for incompressible flow in which the 
assumption of constant properties is valid. For our case/ 
the flow is compressible and the properties/ therefore/ are 
not constant. For applications to Preston-tube data/ the 
properties in Equation should be evaluated at tne wail 
(cone su rr ace) till/ i.e./ 

( 8 . 6 ) 

Tne choice ot wall properties is consistent with Rradshaw 
and Unswor tn*s correlation C201 for comp>r ess ib l e/ turbulent 
houndary layer. 


t D 
w 

— 

4 »w V 


■ F 


AP If 

‘ 4 % V 


8.2 Cnoice of tne Function f 

n atel FI 6 1/1.171 established calibration curves lot the lami- 
nar sublayer/ butter or transition region ana 
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f ully-tuiou 1 Rut. layers. His correlation ier the laminar 


sublayer is 

y* = 0.5 x* + 0.037 (P.7) 

where x ana y are cefined/ lor coupressible flow/ as fol- 
1 ows : 

x * E 1o 9l0 tlHK (B.8a) 

H ww 

y* = lo 9l0 4o ;f ~2 (8.bb) 

p w w 

Alternatively/ x* and y* can be expressed In the lorro 


x* = log 


j * = 


l>£ ■•] 

' S ] 


where Cp - Preston-tuue 

~ ^Pt w ^ ^ ^ Pe 


pressure 


C f = local skin friction cceificient = 


(8.8c) 
(8.8d) 
c cef iicient 


V‘ ,Cj p e U e 


•and 


Pq -= Prestcn-tube Reynolds number 

it was decided to try a straight-line 1 correlation tor 
lamirur ocundary layers in analogy with Equation (8.7)/ 
since it was expected that the behaviour oi a lanirar sub- 
layer Is similar to that oi the laminar boundary layer. 


M.atwt investigation shorten that usin<j a second-order 
function ole not mprove the curve-i ittinq accuracy. 
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8.3 Tho Curve-Fitting Progratu 

The computer program utilized tor the curve fitting is 
called CURF1T. After applying Equation (3.1), the Preston- 
tube pressure data is read in, together with other parame- 
ters like P w / C f / v w , ana U 0 obtained from the Wu ano 

Lock til and the STANS C2J computer programs. 

The probe characteristic lengthy n, has first taken to 
be equal to the height of tho probe. i.e./ 0.0097** (See Fig- 
ure 2). Hut since Fatel*s correlations were based on round 
pcooes/ it was decided to use an equivalent diameter of the 
probe. This was done by assuming the probe face to be an 
ellipse with major ard minor axes 

2a - 0.0138 + 0.004 = 0.0178'* and 
2 D = 0.0097"/ respectively. 

Then the equivalent circle nas an area of 

m n 2 = irab 
4 eq 

Dgq = V4ab' 

D ecf °* 01314" 

The program then calculates x* ard y* fcr each observation 
pointy anc via a curve-fitting package prepared by Or. 
Chandler tl07 at Oklahoma State University/ it fits the 
values of x* and y* to a straight line of the form 

y* = Ax* + K (8.9) 


'Observation points were taxer. 0.5" apart down to the 
end of the laminar portion ol the boundary layer. 


wrier* /> (slope) anc V ( y-intercep t) are constants to be 
determined by the program. 

8.4 Results and Improvements 

Tne resulting straight-line tit to all points was found to 

he 

y* = 0.632 x* + 0.415 (8.10) 

witn a r oot-mean-yguare error 3 of 1.24. This error was con- 
sidered unsatisfactory. 

It was assumed that the reason for this data scatter is the 
correlation model (B.9) does not account for variable pro- 
perty effects. Tnese effects can he accounted lor by intro- 
ducing the reference temperature/ T*/ ir.tc tne correlation. 
At this teirperatur e/ average values for density and viscos- 
ity can ue calculated. Tetervin C183 suggested tnat to 
transfer the incompressible sk in-ir ictior retailor of Ludwig 
and Tillman C19j into Compressible form/ two paraneters neeo 
to be included/ namely M e anc T*/T e . He and numerous other 
investigators have modeled the efiects cf these two parame- 
ters by introducing density and viscosity at a rexerence 
temperature. Although Allen 1.203 selected tne reference 
teirpf ratur e of Sommer and Short (132, we have chosen to use 
chert's formula tor T' as defined m Fguatlon (C.3). Also/ 


•* J e f 1 n e d ay 


E ( y * 5TAN5 ~ y *CURFIT 

»” P° 1nts \ ^STAwT 

No. of points 



tnc USU 


oi the incompressln 1 e RarnoulU 


equation to 


calculate Up^./ Ryuation (8.2), is not accurate. Assuming 


isoeneryetlc flow ^ T t, Pt t,e ^ across tne oounoery Layer. 
up|. can bo calculated more accurately as follows: 

W tl ♦ ) Y/<Y ' 1> 

-2 .2 


or m; 


Pt y-1 


(p pt /p w ) (Y - 1)/Y -i 


( 8 . 11 ) 


and u pt /U e 


_ M pt 

| T pt 

, M pt 

T Pt , 

1 T t,e / 

M e V 

' T e 


T v 

‘t Pt 

t. r 


't,Pt 


**■ V U e 


and x* ts now do fine c as * 
x* 5 log,. (uiLo 2 /4 b 


. Mpt /’ + 111 < \ 1/2 

M e (tT "A m 2 

V + V M Pt, 


t.e 


( 8 . 12 ) 


no • ■ W ' 

Thus, the Improved model is in the form 


( 8 . 13 ) 


y* = A x* + B 1og(T7T e ) + K 


y* = 0.655 x* + 2.095 1 og 1 0 (T'/T e ) - 0.895 


''fie resulting correlation is 

/ ^ I l-l* \ A A A r* 

( 8 . 14 ) 

with an rirs error of 1%. To further ittprove the fitting 
accuracy, a quadratic model or the torn) 


y* = A x* 2 + B x* + C log 1Q (T'/T e ) + K 


was triud. Thu result. Is 

1.2 


y* = 0.273 x* c - 2.618 x* + 1.645 log 10 (T7T e ) + 8.921 (8.15) 

with an rms **rror of 0.054. hquatlon (8.15) can fee written 
in the form 


C f - 6.67 x 10 9 Pw 


10 


,log 2 0 (u Pt°/ 2 V 


0.546 


-5.236 


T e 


(8. 16) 


<>e / U Pt D \ 

\*>7I 

which nas an r ms error of “.85*. Figure 7 shows tne data 
scatter in of . Figure 8 compares the recommended correia- 
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tion (fl.!t>) with the data. The term z* Is defined as 

2 = 0.273 x 2 - 2*618 x ♦ 1.645 log (T*/l ) 

*10 e 

The extraneous data# which appears above the 10% line in 
Figure 1, corresponds to a hach number cf 0.60 and fie ft of 
three and fcur million. It is speculated that these data are 
a result of the iornation of a transonic shock or the stea 
or the flow-angularity probe (e*g. see Reference 8) which 
affects the measured values or P pef and thus P pt . Discarding 
only this particular data/ a new fit results in the 
following equation. 


y* 55 0.0942 x* 2 - 0.438 x* + 2.023 log (T'/T 1+2.272 
10 e' 


( 8 * 1 ?) 

The corresponding rms error in C f is 4.93%. Thus, ‘ Equation 
(8.17) is the recoaaended correlation for relating Cf ano F pi 
within subsonic/ compressible laminar boundary layers. 



Fiqtru ?. Dotation of Predicted Skin-Friction Coelfi 
ciont by Fqn* C^.lb) from Theoretical 
V aiues 






8.5 General Pemarka 


a “ The ^crease o£ U discussed in section 8.3/ resulted In a 
better tit. This can be explained as follows : lr the pro- 
cess of deriving Equation (8.4)/ the non-di mensicna L effec- 
tive center of the probe/ k gf ^ / la assumed to oe unity. 
Patel fl6J and Prozorov C223 and others, e.g. Chue C173, 
have touna that 0 .55<y gff /u<o. o5. Thus, writing Equation 
(8.4) in the lorit 


t-» equivalent to assuming that the average value of y g ^^ = 

1.3 8/2= 0.65 0, or equivalently k _ = 1.3, 

eft 

The better fit is an indication of the strong effect of 
tne prooe geometry expressed by r q . one can also conclude 
tnat k g ^ is a function of Rp. This conclusion was postu- 
lated before by Preston C21J. Patel C163 and Prozorov C223 
shewed that V gff is a function of u pt Q/p . 

b. Although the assumption that K gff is a constant works 

w<i,1/ k eff ia not a COI 'Stant in fact. This can be seen in 
Figure 6, it increases slowly with x. It can be shown that a 
constant k gff requires that the coefficient of x* m the 
correlation be 0.5. The higher coefficient in Equation 
( ;I.10) confirms that k gff is not a constant. Assuoinq a 81a- 
sius type pioiile, it is sho*n in Appendix u that 
ar( x/b) ,337 /( Pa x > - 35S . 

c. Th*» correlation (8.15) u true lor body geometries other 
tnan the cone since it is basec on local variables. It 


account* tot neat transfer conditions since it includes the 
temperature ratio T* /T g • It is also thought to be valid tor 
pressure gradients since it is based on conditions near the 
wall. Thus, It is considered to be a general equation tor 
estimation of skin-friction coefficients in subsonic, lami- 
nar boundary layers* 

<1. Ill an attempt to Improve the curve-tit, ue tried various 
calibration modulo. Among them wore the tolloulng tesults : 
y* ■ 0.102 X* - 0.232 1og 10 M, t 0.815 log, 0 Re ft 

- 0.867 log, 3.458 ( 0 . 18 a) 

with ,jn mis error ot G.27*, 

y* = 0.011 x* - 0.582 log i 0 (1+m/) + 0.481 log 10 Re x 

- 1 .972 log 1Q | - 1.554 (8.18b) 

with an rnss error of 0.06*, and 

log 10 = °- 002 lo 9 10 C(P pt - p 'JqJ + 0.024 log^ (HMJ!) 

+ 0-501 log ]0 (Re ft D £) + 1.688 (8.18c) 

witn an rms error oi 0.02*. Though the accuracy ct the Ht 
became better and better, the dependence on x* anc trie Pre- 
s ton-tube measurements became less and less. This means that 
8TA N r i calculations uete correlated In tnese calibration 
models rather than the experimental data. 

The use of frecstream parameters (M 00 ,Re ft ) ir correla- 
tions («.ifla,b,c) umits thtir use to the 10-cegree cone 
measurements, t.e., the coefficients of N w a no Re ft in these 
correlations are not universal. To correct tor that, the 
local iac.i r.uiubor H g and Re x should be used. 


OUK'.INAU PAOR IS 
O,.' OUAUTY 

ft« Ihe calibration models usoa by Bridshaw anc (Jnsworth/ 

Alleri/ Fenter and Stallmach/ <ina Patel which were reported 

by A lion in his survey report C 20 j were all tried tor the 

present laminar data. It was touno that none ot them was 

competitive with our correlation in terms of the rms error 

In sk 1 n- f r 1 c t l o n coefficient* Allen's 2 nd— degree nodel fit"* 

ted tne laminar data with an rms error in C of 8*64. 

f 

f. Bradshaw and Unsworth C 23 ] have critici2ed Allen's use of 
tiie reference temperature to evaluate density and viscosity 
in the classical law -of -the-wcl 1 • Hather than replace the 

conventional evaluation ot properties at the wall/ we have 
followed the procedure by Tetervin L 183 and others to obtain 
a compressible equation tor C f by simply multiplying an 
incompressible equation for C by the ratio of T*/T raised 
to seme exponent. Here we have determined the exponent via a 
curve fit of the data. Thus/ we are partially accounting tor 
hr adshaw and Unsworth's objection. However. their second 
objection still applies to our analysis in that the refer- 
ence temperature method is based cn zero-pressute-gradlent 
flows and has ar unknown range of validity for flows with 
pressure gradients. 


8.6 Prozorov Correlation 

Assuming a relatively small height of the Preston-tube/ Pio- 
zoruv l 221 expanded the velocity u about the wall using 
Maclaurin's aeries arc reacnec the lollowing simple correla- 
tion for incompressible laminar boundary layers. 


OKI* IiM '•'< 

Ul 40 

(8.19) 

He analytically verified Equation (8.19) for round and rec- 
tanyle openings cl the probe (for which y eff can ue theoret- 
ically calculated). 

Correlation (8.39) has the aovantage that It can be 
uyed foi high pressure gradients*, and the disadvantage that 
Y ff must be known a priori. 

1 1 is also llaiited to incompressible flowa. 

1 r worth irenticniny that Prozorov's paper is the 

only study lound in the literature that ciscusses correlat- 
ing Pteston-tube data with theoretical laminar shear stress. 


" yt 
y off 


r eff 


*A11 the cases investigated in this study tac sirali 
favorable dp/dx. 


CHAPTER IX 


CONCLUSIONS 

!• The Hu and Lock computer program is an accurate and reli- 
able way of calculating tne lnviscid flow field about a 
sharp cone at transonic speeds. 

Irilth the added subroutines, the program is row capable 
oi calculating the inviscid pressure and velocity distribu- 
tion along a conical ray, corresponding to the Preston tube 
survey, tor arbitrary combinations of pitch anc yaw angles. 
Tt also calculates compressible initial protiles based or 
similarity theory anc the supersonic laminar cone rule; this 
information is used to start the boundary layer computa- 
tions. 

The STANb computet code does not work satisfactorily wnen 
the flow angles are significant. tt was found that Its cal- 
culations were insensitive to chances In the flow angles 
when other parameters were kept the same. This limits its 
utility. 

T. It is possible to correlate sxlr. friction and experimen- 
tal Pres ton-tube pressure measurements in the simple form 
( o. 1 7). 

i. The nofi-uimenslonaL effective centur ol the Preston tube, 

* lJ ntt d c ohstant value but rather increases with x 

and decreases with he . 

x 

*il 


CHAPTER X 


supplementary OBSERVATIONS 

1* A 3-diwensional b oundary-layer computer code is needed to 
continue investigation of the role of pitch and yaw angles 
on the correlation of Preston-tube data and skin friction* 

?. The laminar correlation needs to be verified in super- 
sonic flows and also for tree-flight conditions for which 
the wall temperature seldom equals the aaiabatic wall temp- 
erature. 

3. By using the measured Preston-tube pressures at the enc 
ot transition, the correlation oi Allen L203 can be used to 
initiate computation ot the tul ly-deveioped turbulent boun- 
dary layers on the cone. This avoids tackling the develop- 
ment of a skin-tr ic tion correlation tor the bouncary-layer 
transition region until the laminar and turbulent correla- 
tions are establishes 

4. The laminar correlation may be connected with Alien's 
and/or Bradshaw end Unsworth's C20J correlations for turbu- 
lent boundary layers in older to model boundary-layer tran- 
sition. 

r >. In oroet to verify anc make use oi Prohorov's C22J tind- 
itiys, a method is required tnat relates tne Preston-tube 
pressure tc the geottftry of tne proce. Ore way ot coinq this 
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ii; by curve-lit tiny the computed values ot k .. (obtained 

GtT 

troft plots siuilar to Figure o) witi. x and fte • 
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APPENDIX A 


AZIMUTH ANGLE CALCULATION 

In ms appendix are presentee the equations oeveloped by 
Dunn et ai £117 to locate the windward element. Figure 9 is 
a schematic of a typical vehicle at angle ol attack which 
defines the parameters used in this appendix. As illus- 
trated/ the angles of pitch and yaw are. measured with res- 
pect to the freestream velocity vector. It should be noteo 
that the relationship utilized to determine the location of 
the windward element is sensitive to calculation accur_ecy. 
Foe tnis reason/ double precision is used in tne computer 
subprogram ANGLES. The pitch and yaw angles are restricted 
to magnitudes less than 9U degrees. 

Tne first step is to evaluate the angle between the 
vehicle axis ami the resolved yaw vector. This angle will be 
denoted by <£>. 

sin ( <t> ) - c/f 
tan ( <f> ) = c/d 
sin (a) = c/e 
tan ( a ) = c/ a 

sin (a ) / s 1 n ( <M s (c/e)/(c/f) = t/e 
tan i a )/ tan ( <M - (c/a)/(c/d) = d/a 
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c 



Figure y. ycn«matlc cf Flo* Angles 
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cos ( j8 ) = a/d 


ORKiINM. PACI- IS 
OF 1*0015 1)1 ''.M.I'I'Y 


They tan ( 0 ) * (a/d ) tan (o ) = cos ( j8 ) tan (a ) 

tan ( 0 ) = cos ( 0) tan ( a ) (A.l) 

Now the angle a can be calculated as lollaus ; 

ci/f-cos(0), cos (af)= a/f * (a/d)(d/f) 

cos (a ) s cos ( £ ) cos (0)/ <J< cF<9Q (A. 2) 

Kguatiori (A. 2) determines tne angle a which is denoted as 

* 

the effective angle of attack* 

At this point/ he went to find the angle that the windward 
vector makes with tr.c vehicle axis. The 1 ollouiry results 
can be obtained from Figure 9. 


cos t e ) = (0*U')/C0« A) 

0 * A - a sin (o’) 

0*0 = a cos ( ot ) 

tan ( 0 ) = (0*D')/(o»n> 

cos ( e ) = (u*U')/a sin Ca) - (O*o)tan (0)/e sin Ca) 
cos t € ) = a cos (cT) tan { 4> )/a sir (o') 

•*. cos (c)± cot (cT) tan (0) (A. 3) 

an alternate expression for calculating e can he termed by 
substituting Kquatlon (A.l) for tan (0) 

cos (€) = cot («•) COS (0) tan (a) 



APPENDIX H 


CALCULATION OF freestream PROPERTIES 
Values of , q OT anc Re^. are specified for a given uina 
tunnel setting. Prom these values/ all properties of the 
freestream can be calculated as follows : 


t. Uhtain the freestream total pressure P^ ^ as follows : 


q ® = 1 p « u «. 2 = 1 p « M ® 2 (y p ro / P „) 

2 

'CO ' CO ' 


p f" 


* Y M P /2 

l GO co ' 

P 


(B.l) 


t,« 


P co % 


* 


2 q ro 0 


Y-l 2 .Y/(y-1) 0 

+ ia M /Y N 2 


2- O 


CB.2) 

Note that the total pressure ana temperature are constant 
for lsentroplc/ subsonic tlow. 

2. obtain the freestream static temperature as follows : 

2 _ 2 




U p U 

co co oo 


q® 


/ ~yTt~ 




u 2 ^ /TTTT 


2 (T + 198,6) 


2.27 x UT S tJ •* /YTT 

,.(^t.M,r» L R] T2 

\ q„, 2 / m 0 


198.6 = 0 


(R.3) 


Where Sutherland's relation Is used for u • 

”00 
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•1. Using Equations (B.l) and (U.3), p ^ can be obtained from 
the tne perfect qas relation : 

P. " Poo/ R T co IB. 4) 

4. Also the total temperature car be found using the isen- 
tropic r elati on : 

T t = L 0 + ^ mJ) ( B .5 ) 

This procedure is automated in subroutine DtST of the 

exterded Wu and Lock program described in Chapter V. The 
listing or the routine can be found in Appendix F. 


APPKNU I X c 


CALCULATION or INITIAL PROFILES 
Since STANS id a lor ward-marching, l inite-di £1 eter.ce pro- 
yraK, starling profiles of velocity and total enthalpy are 
required to calculate subsequent velocity and total enthalpy 
profiles along the cone* Care shoulo oe taken, therefore, in 
calculating these initial profiles* However, the effect ot 
the starting profile on the calculations oecoroes small after 
a certain developmental distance, as shown m Figure 1* 

The edge velocity distribution can be expressed C13] as 
fol Lows 

u e = c x" <c.n 

Wnere c and n are constants. Fitting equation (C.l) to typi- 
cal edge velocities near the tip ot the cone, as obtained 
from the extended rtu and Lock program 1 , results in 
0 .0047. 

The pressure gradient parameter tor a conical flow, 0 , 
is related to t.he tnviscid velocity distribution 1133 by 

»c ' 3T n - “ 0.003128. 

This in turn corresponds to a wedge flow with 

n W~ V 3 = 0-00157. 


‘Case analyzed here was Case 40.S47. 


Vi 


( )Ul' 

nK I'O 1 " 1 


1'AC.F. I s 

lU’Ml ' 1 V 


‘>3 


Examination of Figure 4-11 of White C133 indicates that the 
solution torf(n) co rresponulng to 0=0 is expected to be 
good. Therefore, the tabulated solution for Blasius flow way 
be used to specify the Initial profiles. The normal dis- 
tance y c can be calculated now trom 

y c * 

^c*v[r^r ] 5 < c - 2 > 

Where v' is the kinematic viscosity evaluated at the refer- 
ence temperature T' as will be shown now. One can obtain 
Equation (C.2) using Mangier transformation.— 


An expression for the reference temperature across the 
boundary layer is given by Eckert's formula C133 : 



t'= T(0 . 5 + 0.039 M 2 + 0.5 T /T) 

w 

(C* 3 ) 

Where 

T = V(1 *!^!m 2 ) 

(C.4) 

and 

T w ’ T aw ’ ^ + r M 2 ) 

(C.5) 

Where r 

1 /2 

- (Pr) for laminar boundary layer, and 



1 /3 

= (Pr) ' for turbulent boundary layer. 



The values of local Mach rumbers and the total tempera- 
ture Tj. are calculated by the extended Wu and Lock program. 
Prendtl number is taken to be 0.72 for air. 

Now At can te celculated using Sutherland's relation 


p = 2.27 x 10~ 8 


(T) 1 ,a 
r + 198.6 


(C.6) 


lo obtain v' , p' Is calculated using the perfect gas 

relation 


V R 


T' 


IC.7> 


( > <s 

tfH ' at,ltic Measure t' w la calculated ir. tac main Wu 
«nc1 Lock program. 

Therefore, from (C.6) and (C.7), y' can be calculated/ 

v ' = vh' (C.R) 

Substitution in {C.2), yields a table of y vs. 7 j D . Prow the 

niaslus solution tor f = u/« e , „e can obtain a table of u 
va. tj b . Thus tne initial velocity profile is specified. 

Thu total enthalpy Is defined as 

2 

h t = h + r!Hr 


' c p T + 2 9, j 


for Cp = constant (C.9J 


n.e distribution ot T through the boundary layer may be 
approximately expressed C13j as 

T * T + (T - T ) -*d — r u 2 

w v aw V U e 2 c p g c J 

Substitut.ion into (C.9) gives 

h t = c [T = (T - T ) -'Ll + (l - r) u 2 

t pw v aw w ; U J 2 q J 

e 

tfle df >su»pt ton 1 = T , this equation reduces to 

vv aW 


h = c T + i l — : IjtL 
t C P 'w 2 g c J 


(C.10) 


T w calculated via the Wu aria Lock program, pr.d tne ini 
tiai stagnation proliie was defined by Equation (C.lO). 


APPENDIX 0 


FUNCTIONAL DEPENDENCE CF THE EFFECTIVE 
CENTER OF THE PROBE 


For simplicity/ i»e hill derive an expression tor for 

incompressible flow over a flat plate. The correlation 
(8.14) reduces lr this case to 


r r 0.655 D 0.69 
L f~~ L p K D 


(D.l ) 


-.'here C f = t w /( 0.5 pU , 

C p = AP / ( 0. S pU 2 )/ 

^ = U U/v s <U x/v ) n/X s Re x D/x. 


Since AP = ~ P u p t 2 


C can be written In the torm 

r 

c p s (u pt /u > 2 ■ <‘ , J 2 atn 


eff 


wnere £• - 1st derivative ot the Blasius function w.r.t, 
and 


h/ 


n eff = Vf H® 72 v x ~ k eff 0 /Ri 7 /X - 
"incc the height of the probe is very small (0.0GS7**)/ all 
the laminar boundary-layer uata was obtained hlthm the 
lower oi the layer thickness. In this region f»~T? is 


valid. 

Therefore/ Cp 
Substl tu tiny 


can be expressed as 

C p~ n eff u fc2 eff ° 2 Re x /x2 * 
relation (D.2) into (D.l) give; 


r 1,1 .31 Dfi -0.035 /Xx -0.441 
l f~eff Re x ( D' 


( 0 . 2 ) 


(D.3) 


The well-known relation for C in this case Is 

f 

C f ~Re x -0,5 (D.4) 

Comparing (D.3) and (b.4) y the ioi.' wing equation is 
obtained tor k e ff* 


k n -0.355 ,XvO. 337 
k eff Re x 


(D .5a) 


Or alternatively/ Re y " 0 ’ 018 R r "°‘ 337 


(0.5b) 


'eff “x ’'D 

Again/ relations (D.5a/b) are only valid for incompressible 

flow over a flat plate and ctre presented nere only to 

demonstrate that y fi ff is not a constant. 
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RAW DATA USED FOK SKIN-FRICTION 
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